Spontaneous mass current and textures of p-wave superfluids of trapped Fermionic 

atom gases at rest and under rotation 
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It is found theoretically based on the Ginzburg-Landau framework that p-wave superfluids of 
neutral atom gases in three dimension harmonic traps exhibit spontaneous mass current at rest, 
whose direction depends on trap geometry. Under rotation various types of the order parameter 
textures are stabilized, including Mermin-Ho and Anderson- Toulouse-Chechetkin vortices. In a 
cigar shape trap spontaneous current flows longitudial to the rotation axis and thus perpendicular 
to the ordinary rotational current. These features, spontaneous mass current at rest and texture 
formation, can be used as diagnoses for p-wave superfluidity. 

PACS numbers: 03.75.Ss, 67.85.-d 



Superfluids with multi-component order parameter 
(OP) form a unique physics branch, including from 
spinor BECP, 0, U 0] in cold Bosonic atom gases with 
^^Na or ^^Rb, strong interacting Fermionic liquid '^He 
atoms S 0, S|> some of the heavy Fermion supercon- 
ductors, such as UPtaj^, It is further extending to 
color superconductivity in dense quark-gluon plasmas in 
high energy physics [III [l^. This is unique because one 
can expect rich topological defect structures or vortices 
and explore a new phase of matter interesting from a fun- 
damental physics view point. Some of vortices, which are 
able to accommodate Majorana zero mode at a core may 
be useful in quantum computing [13]. 

Recently there have been much attention, in this 
respect, on p-wave resonance superfluidity made by 
Fermionic alkaline atom gases, such as ^Lilll [H and 
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40k[T1, 

[20j . Since experiments for achieving 



both experimentally and theoretically [lil 

wave resonance 
superfluidity are steadily progressing [21|, it might be a 
good timing and also necessary to consider the generic 
properties of p-wave superfluidity both at rest and under 
rotation in order to help detecting its superfluidity where 
non-trivial topological structures, or vortices are a hall- 
mark to identify it. Namely, the textural structure of the 
three component order parameters spanned by three ba- 
sis functions Px, Py and Pz- The order parameter space 
consisting of three components is analogous to that of 
the superfluid '^He, in particular, the A phase[^ H, 0| 
where the OP is described by a tensor A^.i = df^Ai (/z, 
i = x,y,z). The d-vector (Ai) describes the spin (orbital) 
state of a Cooper pair. Since in the p-wave Feshbach res- 
onance superfluid under a magnetic sweeping the spin 
degrees of freedom is frozen, only the orbital degrees of 
freedom is active, namely the OP is characterized only 
by Ai. In a sense our p-wave superfluid is analogous to 
the "spinless" superfluid ^He-A phase 22 1. 

A critical difference between superfluid '^He and p-wave 
resonance superfluid of atom gases lies in the boundary 
condition: Atom gases are confined three-dimensionally 
(3D) by a harmonic trap, which is easily controlled, re- 



sulting in various shapes, such as cigar or pancake shapes. 
As we will see soon, the trapping potential gives rise 
to an important handle to control the the 3D texture 
structure. Indeed, the 3D trapping structure constrains 
possible textures as a whole. This feature is usually 
absent in superfluid '^He where 2D textures, such as 
Mermin-Ho [23 or Anderson- Toulouse-Chechetkin^2j| are 
discussed. They are essentially 2D textures. Here we are 
interested in finding a truly 3D object. 

The dipole-dipole interaction between two alkaline 
atoms acts to split the relative orbital state for two par- 
ticles, depending upon the projections of the orbital an- 
gular momentum, m; = ±1 and mi = 0. This results 
in breaking the degeneracy among A^ ± iAy and Az, 
or px ± ipy and Pz are non-generate. This splitting is 
estimated to be large for ""^K by Cheng and YipLlQ.] evi- 
denced by the clear difference in the Fershbach resonance 
magnetic fields (the splitting field=0.47±0.08G)[l3. For 
^Li it could be small where there is no clear observation of 
the resonance splitting at a magnetic field H = 158.5(7)G 
in Univ. Tokyo experiment so far 21| . In this paper we 
take its splitting as a parameter to extract the generic 
features of the p-wave superfluids. Our purposes are to 
find the possible 3D texture formed by the /-vector de- 
fined shortly in 3D harmonic trap and help identifying 
the p-wave superfiuidity. 

The superfluid condensate is described as A = 
AxPx + AyPy + AzPz = A+P+ + A-p- + Aopz with 
P± = TiPx ± iPy)/V2, A± = TiAx T iAy)/V2, 
Aq — Az- A± and Aq are the three OP's for the p- 
wave superfluid. The Z- vector is deflned as Ix{y) = 
V2Re{iA+ + A-)A*o}/\A\^ly{r) = V2Im{{^A+ + 
A_)A*o}/\A\Mz{r) - (|A+|2-|A_|2)/|A|2. Thisf-vector 
now fully characterizes the p-wave superfluid and its spa- 
tial variation is called the Z-vector texture or simply tex- 
ture. 

Here we employ the Ginzburg-Landau (GL) frame- 
work, which is free from the microscopic details and 
which is able to lead to generic topological structure of 
the problem. Namely, we start with the GL free energy 
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functional which is expanded in terms of the OP's A±{r) 
and Ao{r) up to the fourth oder, 

f fgrad ^ fbulk ^" fharvaonic + fc 
fgrad = i^l (9* A* ) (9. A, (9* A* ) (9, A.)+7^3 (5* A* ) (3, 

fbulk = -ao(l - U)A*A, + P24A*A,A*Aj+P3A*A*AjAj 

f harmonic = ]^muj']_{p^ + \^ Z^)\Ai\^ 

where ti — T /Td (Td is the transition temperature for 
i-component) , di = Vi — x r)i, and = + . 
We consider £7 || z. The anisotropy of the harmonic 
trap is expressed as A = cj^/wx. The GL parameters 
ao, P2i = P2 + /34, and Ki = K2 = K3 = K are 
taken as those estimated by the weak couphng approx- 
imation, assuming the Fermi sphere^: ao = 7V(0)/3, 
= P3 = Pi = 7C(3)iV(0)/120(7rfcB7;)2 and K = 
7C{3)N{0){hvF)^ /2A0{TrkBTcy where N{0) is the density 
of states at the Fermi level, and vp is the Fermi veloc- 
ity. The weak coupling approximation should be a good 
guide to understand the generic properties of the p-wave 
superfluids of atom gases because it has applied success- 
fully even to liquid ^He of strong interacting Fermions 
with only additional small strong corrections [sj . 

It is interesting to notice that the centrifugal poten- 
tial leads to non-trivial form fc = — ^m^l^ {\Ap\^ + 
3|Aep| + A^p) with Ap = A^cos9 + Aj^sin6', Ag = 
—Ax sin 9 + Ay cos 9 because the OP label implies the or- 
bital angular momentum, a feature absent in the spinor 
BEG. The extra factor 3 in the above becomes impor- 
tant when evaluating the critical angular velocity flcr 
above which the superfluid turns to be normal. That 
is, flcr — uj±/V^ is greatly reduced from the usual case 
(flcr = L0±). The trapping potential f harmonic acts as 
lowering the transition temperatures. As mentioned, the 
dipole interaction causes the splitting of the transition 
temperatures into the two groups Tex = Tcy and Tcz- 
We introduce a = Tex /Tcz (0 < a < 1) which indicates 
the degrees of the broken symmetry of the system, that 
is, the three components are completely degenerate for 
a = 1. When a ^ 0, the one component Ao{^ 0) or 
scaler superfluid tends to be realized. 

Before going into the confined system, we first con- 
sider an infinite system. The phase diagram shown in 
Fig.l consists of the three phases: The lower A phase is 
described by a chiral OP Aq + iAj (j = x or y) where 
the time reversal symmetry is broken. The high tem- 
perature B phase is described by Aq which is the one 
component scaler superfluid. In the following we exam- 
ine the A phase in the confined geometries. 
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FIG. 1: (Color online) Phase diagram of the chiral p-wave 
pairing state in infinite system in temperature {tx ~ ty) ver- 
sus the anisotropy a = T^x/Tcz- N: normal state, B: single 
component state Aq and B: the chiral state Aq + iAj {j — x 
or y). 



We have found stationary solutions by numerically 
solving the variational equations: bFjbAiiv) — {i — 
±, 0) in three dimensions. For a cigar shape trap, 
80 X 80 X 120 meshes are taken with the cloud sizes 
Rx = Ry = 10^0 and Rz = 50^o- The coherent 
length ^0 = -64:9 hvp/2TTkTc. For a pancake shape trap, 
100 X 100 X 80 meshes are taken with the cloud sizes 
Rx — Ry = 30^0 and Rz = lO^o- We started various ini- 
tial configurations, which are uniform, include singular 
or non-singular vortex, to make sure that the resulting 
texture is energy-minimum. Throughout the paper the 
temperature tx — ty = 0.4 and a = 0.9 are fixed for the 
cigar trap (A = 0.2). For the pancake trap (A = 3.0) 
tx — ty — 0.5 and a — 0.95 are fixed 

It is noted that in our problem the boundary condition 
is essential in determining a stable texture. In superfiuid 
■^He-A phase the ^-vector is always perpendicular to the 
hard wall so that the perpendicular particle motion is 
suppressed. In other words, the point nodes situated to 
the ^-vector direction touch the hard wall so as to min- 
imize the condensation loss at the boundary^. Now in 
our harmonic trap where the condensation density gradu- 
ally decreases towards the outer region, the Z-vector tends 
to align parallel to the circumference. This orientation is 
advantageous because the condensation energy is maxi- 
mally gained by letting the point nodes move out from 
the system. We can expect quite different situations for 
cigar and pancake shapes, which indeed leads to quite 
distinctive 3D textures. 

We start out with the stable texture at rest for a cigar 
shape trap with A = 0.2. Atfi^O the ^-vector is shown 
in Fig. 2. Figure 2(a) displays the amplitude distribution 
of the ^-vector. It is seen that the amplitude |Z| is max- 
imal in the central region and towards the outer regions 
1^1 becomes smaller gradually. At the top and bottom 
ends the polar state is realized where the ^-vector van- 
ishes. The three cross sections are shown in Figs. 2 (b), 
(c) and (d). In Fig. 2(c) which corresponds to the middle 
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FIG. 2: (Color online) Stable texture at rest for the cigar 
trap A = 0.2. (a) Distribution of |/| in z-x plane, (b)-(d) 
Three cross sections shown in (a) for Iz (color) and Ix and 
ly components (arrows), (e) Spontaneous circulating current 
flows in z-y plane along the z direction. 
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cross section the /-vectors lie in the x-y plane, showing 
a stream line type pattern in which the /-vectors follow 
the circumference as if the water streams along the circu- 
lar boundary. In the outside of the condensate, a unseen 
sink and source of the /-vector exist so that two imagi- 
native focal points situated outside. The left (right) one 
corresponds to a source (sink) where the /-vector appears 
(disappears). This stream line like texture is contrasted 
with the so-called Pan- Am type texture in superfluid '^He 
A phase [6] where the /-vectors tends to point perpendic- 
ular to the wall due to the boundary condition. In the 
upper (Fig. 2(b)) and lower (Fig. 2(d)) cross sections 
the stream line type textures maintain, but Iz compo- 
nent appears additionally. 

The associated mass current structure is depicted in 
Fig. 2(e). The jz component shows the circulation current 
along the z axis. Since il = 0, this circulating mass cur- 
rent is spontaneously generated. This non-trivial mass 
flow is explained in terms of the so-called bending cur- 
rent: cx V X /. This is not usual direct current due 
to the OP phase modulation because of = 0. There 
is no phase twisting here. This is due to the /-vector 
bending. In the central cross section in Fig. 2(c) the /- 
vector in-plane bending (V x l)z produces the perpen- 
dicular current jz- However, at the upper (lower) plane 
Fig. 2(b) (Fig. 2(d)) the mass current acquires the jx and 
jy components because of the non-vanishing Iz compo- 
nent appearing there. Therefore, the perpendicular cur- 
rent at the center bends around so that the total mass is 
conserved as required. It is clear from Fig. 2(e) that the 
mass current circulates perpetually along the z direction 
parallel to the long axis of the trap. This result is non- 
trivial and remarkable to manifest itself the topological 
nature of the / texture. 



FIG. 3: (Color online) Texture change with rotation speeds 
Q, relative to the in-plane trap frequency Central cross 
sections are displyed. 1^ (color) and Ix and ly components 
(arrows). SI = corresponds to (c) in Fig.l. As increaing Q 
/-vectors acquire the z component (color changes). At 51 = 
0.4:Lo± two Mermin-Ho vortices enter from y direction seen as 
yellow objects. At f2 = 0.5uj± four MH's are visible and the 
condensate expands and distorts. 

Under rotation, the above 3D texture deforms contin- 
uously and smoothly so as to accommodate the in-plane 
circular direct current of the usual type. It is seen from 
Fig. 3 where the central cross sections of the cigar under 
rotation are displayed that as rotation increases, (1) the 
/-vectors completely in the x-y plane pointing x direction 
acquire the z component as seen by color change from 
green to blue. (2) Above a certain rotation (fi = 0.4:UJ±) 
a pair of Mermin-Ho (MH) like vortex[23 enters from 
the y direction, where at the core the /-vector point- 
ing to z direction flares out circularly and at far sites 
it lies almost in x-y plane. This vortex is nothing but 
Mermin-Ho structure. (3) Upon further increasing rota- 
tion (n = Q.5u!±) two pairs of the MH vortex appears. 
(4) Gradually and concurrently the background /-vectors 
point to the negative z direction. (5) Moreover, the con- 
densate profile itself deforms and deviates from circular. 
Simultaneously it spreads out due to higher rotation, seen 
as extended area for H. = 0.5a;^ case. 

Figure 4 shows a different view of Fig. 3 where the z-x 
cross-section is displayed. At $7 = it is clearly seen that 
all /-vectors point to the x direction. As fl increases, the 
down-ward Iz component appears, which is responsible 
for the counter-clock- wise circular mass current. At the 
top and bottom ends the polar core vortices appear as 
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FIG. 4: (Color online) Corresponding z-x cross sections to 
Fig. 3. The Z-vectors are almost in x-y plane at = 0. Under 
rotation two polar core vortices shown in stars appear. At 
il = O.StJx two of the Mermin-Ho vortices out of four can be 
seen as indicated by red arrows. 
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FIG. 5: (Color online) Stable texture in the pancake trap with 
A = 3.0 at rest, (a) The Z-vector pattern in z-r plane. It is 
axis-symmetric around z. (b) Spontaneous mass current in 
x-y plane at z = 



Z-vector decreases towards the outside. At the center the 
Z-vector is completely locked to the negative z direction 
and towards the outside it becomes twisting, which ulti- 
mately generates the bending mass current in x-y plane 
as shown in Fig. 5(b). However, this in-plane circular 
current diminishes further outside because the absolute 
magnitude of the condensate decreases. Therefore, the 
current maximum occurs in the intermediate circular re- 
gion. 

Under rotation, in stead of the MH vortex in the cigar 
case, a pair of the Anderson- Toulouse-Chechetkin(ATC) 
vortex[24| appears one by one as increasing rotation. 
This is understood because the overall Z-vectors tend to 
align to the negative z direction from the outset, thus the 
vector of the core with the positive z direction changes 
over completely the surrounding /-vector with the nega- 
tive 2 direction. This object is nothing but ATC vortex. 

It might be interesting to point out that the direction 
of the spontaneously generated mass current at rest is 
always perpendicular to the direction of the majority l- 
vectors, that is, in the cigar (pancake) case the /-vectors 
point to x-y (z), so the mass current flows to the z (x-y) 
direction. This implies that in our system the shape is 
decisive in understanding and controlling the physics of 
the textures. 

In conclusion, we have found stable /-vector textures 
at rest and under rotation for p-wave chiral superfluid 
to be realized by using Feshbach resonances of neutral 
atom gases. These textures are quite generic independent 
of the anisotropy a, provided that temperature is low 
enough, or the system is in the A region in Fig.l. The 
spontaneous mass current and the /-vector textures can 
be used as diagnoses to detect and characterize p-wave 
superfluidity. 

We thank T. Ohmi, M. Ichioka, T. Mizushima, and T. 
Kawakami for useful discussions. 



indicated by star marks. The = 0.5tJj^ case shows the 
side view of the MH vortex: The red lines indicated by 
arrows correspond to the cores of MH, which bend out- 
ward, inducing the profile modification of the condensate 
near the arrows. 

Let us now discuss the pancake shape trap. From Fig. 5 
where the resulting /-vector texture (Fig. 5(a)) and asso- 
ciated in-plane mass current (Fig. 5(b)) at = are dis- 
played, it is seen that most /-vectors point to the negative 
z direction except that near the upper and lower surface 
region /-vectors acquire x-y component. The boundary 
condition enforces the /-vector parallel to the surface. It 
is especially strong when the curvature of the surface is 
large. Thus, in this case (see Fig. 5(a)) the left and right 
ends force the /-vectors point to z direction, determining 
the overall /-vector configuration to the negative z direc- 
tion, even in deep inside vectors. The amplitude of the 
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